SHELL MODEL APPLICATIONS TO N = Z NUCLEI* 

Gabriel Martinez-Pinedo 

Department fiir Physik und Astronomic, Universitat Basel, CH-4056 Basel, 

Switzerland 

(Received April 27, 2001) 

This manuscript reviews recent results of large scale shell model cal- 
culations O'o. N — Z nuclei. The topics discussed include the role proton- 
neutron pairing on the binding energies of p/-shell nuclei and his influence 
on backbending rotors; coulomb energy differences and their relationship 
with alignment. Proton-neutron correlations are found responsible for the 
characteristic features of the isotopic shifts in calcium, the parabolic de- 
pendence on A and the prominent odd-even staggering. 

1. Introduction 

In the recent years, thanks to the enormous progress done in computer 
capabilities and programming techniques, shell-model calculations, based 
on microscopic effective interactions and in model spaces that include all 
the basic physics ingredients, have become available. The spherical shell 
model approach has provided an excellent description of many properties, 
from level schemes to rare processes as double beta decays, and has given 
a unified description of the single particle and collective degrees of freedom 
of the nucleus. Taking advantage of the availability of detailed microscopic 
calculations we will study the role of proton-neutron correlations in = Z 
nuclei. We will discuss the effects of pn pairing on the binding energies of 
p/-shell nuclei and his influence on the rotational properties taking ^'^Cr as 
an example. Next we will examine how the alignment of particles evolves 
as a function of angular momentum and his relationship with the observed 
coulomb energy differences in mirror nuclei. Finally, we will show how cross 
shell proton-neutron correlations are responsible for the isotope shifts in 
calcium. 



* Invited talk presented at the "High Spin Physics 2001" NATO Advanced Research 
Workshop, dedicated to the memory of Zdzislaw Szymanski, Warsaw, Poland, Febru- 
ary 6-10, 2001. The results presented in this talk have been obtained in collaboration 
with E. Caurier, K. Langanke, F. Nowacki, A. Poves, J. Sanchez-Solano, P. Vogel, 
and A. Zuker. 
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2. Pairing and N = Z nuclei 

The study of the isovector pairing interaction among hke particles is one 
of the classical themes of nuclear physics. Proton-neutron pairing has been 
much less studied, in particular its isoscalar part. The advent of high-spin 
physics and unproven spectroscopic data in heavy N ^ Z nuclei has again 
brought an old question: how does pairing evolve with angular momentum? 
We shall examine these issues in what follows. To start we select as our 
r = and T = 1 pairing Hamiltonians those extracted in [1] from the 
realistic G-matrices. We keep their notation and call POl and PIO the 
isovector and isoscalar L = pairing Hamiltonians. In the absence of any 
other interaction their spectra are given by: 

„\in-Vs){^n + 6-n-Vs) , t{t + 1) T{T + 1) 
Epoi = -G + 



^\{n-vt){An + Q-n-vt) , s{s + l) 5(5 + 1)" 

^Pio = -G y + 

where G is a coupling constant, $7 is the maximum number of L = states 
in the valence space, n the number of valence particles, Vs and vt the singlet 
and triplet seniorities, t and T the reduced and total isospin, and s and S 
the reduced and total spin. It is interesting to note that while the isovector 
pairing favors states with good isospin the isoscalar pairing favors states 
with good spin. Due to the presence of a strong spin orbit term in the 
nuclear force the presence of a L = isoscalar condensate in the ground 
state oi N = Z nuclei is very unlikely [2] . 

In order to determine the role of the pairing operators in the behavior of 
different nuclei and different physical quantities, we make first a reference 
calculation using the interaction KB3 [3] and then we subtract from KB3 the 
isovector or isoscalar pairing Hamiltonians and make the calculations with 
the new interactions KB3-P01 and KB3-P10. We obtain the effect of each 
pairing channel by direct comparison with the reference calculation. The 
value of the coupling constant G is obtained from the numbers in Table 1 
of ref. [1]. We use G = -0.295 MeV for POl and G = -0.459 MeV for PIO. 

Figure 1 shows the contributions to the ground state binding energy 
from the isovector pairing (labeled POl) and from the isoscalar pairing (la- 
beled PIO) for several isotopic chains. The figure shows the strong odd-even 
staggering of the POl points, mostly suppressed in the PIO. The only little 
surprise is that for A=46 and 50 moving from T = 1 to T = 0, not only 
the contribution of POl decreases, but also the contribution of PIO. This 
may explain why the ground state of ^^V and ^°Mn has T = 1 instead of 
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Fig. 1. Pairing contributions to the ground state energies of the A = 46, 48, and 50 
isobaric chains as a function of the isospin (T) 



T = 0. The monopole part (symmetry energy) part of KB3 will put the 
centroid of the T = states lower than the centroid of the T = 1 states by 
about 1.3 MeV; on the other hand, the total pairing contribution to T = 1 
is larger by nearly 2 MeV than the contribution to T = 0. Therefore, it 
is the T = 1 state that becomes the ground state of the odd-odd = 
nuclei. 

To study the evolution of the pair content of a nucleus as the rotational 
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The effect of the isoscalar and isovector pairing in the gamma energies 
yrast hne in '*^Cr. b) Gamma energies obtained with the Gogny force in 
Model and cranked Hartree-Fock-Bogoliubov approaches. 



frequency increases we have chosen °Cr, which is the most representative 
example. In figure 2a we present as a backbending plot the effect of the 
pairing correlations, basically a change in the static moment of inertia. For 
states below the backbending - i.e. for those that can be viewed as proceed- 
ing of the same intrinsic, well deformed state - the wavefunctions with or 
without pairing have the same structure. Their i3(£'2)'s and quadrupole 
moments are equal within a few per cent and their overlaps are always 
better than 95%. This means that pairing does not affect the quadrupole 
properties that have somehow reached saturation in the deformed regime. 
This suggest that the discrepancies found in the comparison of the CHFB 
and SM results in ref. [4] are due to deficiencies in the treatment of pairing 
in the mean field description. To clarify this point we have used the ma- 
trix elements of the Gogny force [5] (obtained [6] using the wave functions 
computed in a spherical Hartree-Fock calculation with the same force) . The 
results are plotted in figure 2b. It is clearly seem that the shell model cal- 
culation with the Gogny force now gives a much better moment of inertia 
than the CHFB calculation with the same force. 



3. Coulomb displacement energies and alignment 

In the previous section we have study the pairing correlations as a func- 
tions of the rotational frequency. However, traditionally phenomena as the 
backbending of the appearance of yrast traps are related to the alignment 
of particles. In ref. [7] it was show that the 12"*" yrast trap in ^^Fe can be 
related to the alignment of two pair of particles. We can use a schematic 
Hamiltonian to "count" the number of aligned particles: 




i^align = E ^IsiJT) ■ Zrs (JT) , (3) 

rs 

where the sum runs over different orbits (rs) and Zj.g{JT) is the creation 
operator for a normahzed pair of particles in the orbits rs with angular 
momentum J and isospin T. As the conmutator [Zj.g{JT), Zrs{JT)] is not a 
c-number the operator above is not a true number operator however it can 
provide qualitative information about the change in the number of pairs. 
In figure 3 we present the expected value of -ffaiign ( J = 6 is the maximum 
spin with T = 1 in the pf shell) for the yrast band of ^^Fe. We include also 
the expected value of POl using G = 1. For low spin states (J = 0-4) ^'^Fe 
presents a rotational spectra. As can be seem in the figure the rotations is 
collective, no alignment of particles takes place (flat -ff align)- For J = 6 a 
pair of particles is aligned. As J = 6 is the maximum possible spin with 
only a J = pair broken a second pair is broken producing a sudden drop 
in POl. A progressive increase in -fTaiign takes now place and for J = 12 
two pairs of particles are fully aligned. Note that this sudden alignment 
of particles could be responsible for the very small B^EA) (5 x 10~^ W.u.) 
recently measured in a study of the decay of the 12+ state [8]. For spins 
higher than 12"*" there is a clear correlation between the points of maximum 
alignment (in the p/-shell there could be only 2 pairs with J=6 and T=l) 
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and the breaking of a J = pair. 

Gamma-ray spectroscopy of liigh spin states in mirror nuclei has been 
revitahzed in recent years through the development of large gamma-ray 
spectrometers. Such spectrometers have yielded a wealth of new informa- 
tion on high spin states in N = Z nuclei and proton-rich nuclei in the 
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shell. The small differences (normally a few tens of keV) observed in the ex- 
citations energies are mainly due to the isospin symmetry breaking Coulomb 
interaction. In recent works [9, 10], the high spin states of T = 1/2 mirror 
nuclei ^''Cr/^'^V and ^^Mn/^^Cr have been studied. It was found that the 
CED's are extremely sensitive to the structure of the nuclear wave functions. 
In recent experiments [11, 12] high spin states up to the band termination 
in the mirror pair ^-'^Fc/^-'^Mn have been observed. The measured CED's are 
shown on figure 4 where they are compared with shell-model calculations 
that include the Coulomb matrix elements denoted "A42" on ref. [9], they 
were obtained from the A = 42 mirror pair. It is seen that there is an 
abrupt change in the CED's at J = 17/2~. The effect is also present in the 
calculation, that shows exactly the same trends as the experiment, although 
with enhanced values. (Sec ref. [13] for a quantitative reproduction of the 
measured CED's.) The large increase in the CED can be interpreted as 
due to the alignment of a pair of protons in ^^Fe which reduces the overlap 
of their spatial distributions and causes a corresponding reduction in the 
Coulomb energy. In ^^Mn, the odd proton blocks this alignment and a pair 
of neutrons is aligned, with no resulting coulomb effect. Beyond J = 17/2~, 
the protons start aligning also in ^-"^Mn and therefore the CED's approach 
zero at the band termination. In order to make more visible we have com- 
pute the expected value of the proton contribution to /^aligm defined in 
equation (3). Figure 4c shows the difference between the expected values 
for states of ^^Mn and ^^Fe. 

It is interesting to note the similarities between the alignment observed 
in ^^Fe and the mirror pair ^^Fe/^^Mn. In fact the aligned J = 17/2~ state 
in ^^Fe/^^Mn could be seen as resulting of the couphng of a j = 7/2~ hole 
to the yrast trap 12+ state in ^^Fc. Similarly to what happens in ^^Fe where 
a very small B{E4) decay probability has been measured, the 17/2~ state 
in ^^Fe/^^Mn has also very small electromagnetic transitions [14]. 

4. Cross-shell correlations: isotope shifts in Calcium 

The appearance of shell gaps associated with magic nucleon numbers 
is one of the cornerstones of nuclear structure. However, it has become 
increasingly evident in recent years that these magic numbers, and the cor- 
responding shell closures, might get eroded with increasing neutron excess. 
A prominent example is the magic neutron number N = 20 which vanishes 
in proton-deficient nuclei with Z < 12. This erosion of the shell closure has 
been related to cross-shell proton-neutron interaction which correlates the 
2si/2 and orbitals with the I/7/2 and 2^3/2 orbitals [15]. Similar cross 
shell correlations are also responsible for the superdeformed band seen in 
^^Ar [16] (see also C.E. Svensson contribution). If cross-shell correlations 
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are indeed the mechanism for the shell erosion, then first indications are al- 
ready visible in the stable calcium isotopes. Our argument is based on the 
understanding and explanation of nuclear charge radii, (r^), in the calcium 
isotopes. The isotopic shifts show a characteristic parabolic shape with a 
pronounced odd-even staggering when neutrons fill the /7/2 orbit and the 
mass number changes from A = 40 to A = 48. 
The nuclear charge radii can be written as: 

= rc(mean field) -|- (correlations). (4) 

Mean field calculations, which usually aim at describing nuclear masses, 
deformation parameters, and radii over a large region of nuclear masses and 
charges, cannot account for the details of the calcium isotope shifts. The 
dependence of (r^) on A is usually featureless. Some of the approaches, 
however, arc able to account at least for the near equality of (r^) in ^'^Ca 
and ^^Ca. To this category belong calculations based on the Hartrec-Fock 
method with Skyrme interactions [17], relativistic mean field methods [18], 
and the extended Thomas-Fermi model with the Strutinsky-integral [19]. 
Then we will assume that r^(mcan field) in equation (4) is constant for 
all the calcium isotopes and determine the contribution of the correlations 
using shell-model calculations. 

Since here we are interested in the description of calcium isotopes, it 
is imperative to include states in the vicinity of the N = Z = 20 shell 
boundary. Therefore, the chosen valence space consists of the ^3/2; Si/2) 
fj/2, and P3/2 subshells for both protons and neutrons. (Thus ^^Si represents 
the inert core.) Our calculations [20] reproduce the energies of the intriidcr 
states in Sc and Ca, as well as the energies of the low-lying 2+ and 3~ states 
(and B{E2)) in the even Ca isotopes. 

Due to the configuration mixing across the Z = 20 shell boundary, 
protons are lifted from the sd to the fp shell, resulting in a increase of 
that assuming harmonic oscillator wave functions is given by: 

Sr'c = ^n%, (5) 

where Z = 20, h is the oscillator parameter which we assume constant for 
A = 40-48. (Notice that this assumption is supported by the fact that 
is constant in mean field calculations.) n^^ is the number of protons in the 
fp shell that we obtain from our shell-model calculations. 

Figure 5 compares the computed isotopic shifts with data. The trends, 
parabolic shape and odd-even staggering, are clearly reproduced but the 
magnitude of the calculated shifts is smaller than the experiment suggests. 
That could be due to some small dependence of b with A or the fact that 
we have neglect the (^5/2 and /5/2 orbitals in our model space. However, 
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Fig. 5. Isotope shifts in calcium. The experimental data (circles connected by a 
solid line) and the shell- model results (stars connected by a dashed line) are shown. 



our results clearly show that cross shell correlations are responsible for the 
observed isotopic shifts in the calcium isotopes. 
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